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Abstract 

The Higgs bosons of the MSSM can be produced from the decay of SUSY particles. We 
review the evalulation of two decay modes in the MSSM with complex parameters (cMSSM). 
The first type is the decay of the heavy scalar top quark to a lighter scalar quark and a Higgs 
boson. The second type is the decay of the heavy chargino to a lighter chargino/neutralino 
and a Higgs boson. The evaluation is based on a full one-loop calculation including hard 
QED and QCD radiation. We find sizable contributions to many partial decay widths and 
branching ratios. They are roughly of 0(10%) of the tree-level results, but can go up to 
30% or higher. These contributions arc important for the correct interpretation of scalar 
top quark decays at a future linear e+e~ collider. 
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The Higgs bosons of the MSSM can be produced from the decay of SUSY particles. 
We review the evalulation of two decay modes in the MSSM with complex parameters 
(cMSSM). The first type is the decay of the heavy scalar top quark to a lighter scalar 
quark and a Higgs boson. The second type is the decay of the heavy chargino to a 
lighter chargino/neutralino and a Higgs boson. The evaluation is based on a full one- 
loop calculation including hard QED and QCD radiation. We find sizable contributions 
to many partial decay widths and branching ratios. They are roughly of 0(10%) of the 
tree-level results, but can go up to 30% or higher. These contributions are important 
for the correct interpretation of scalar top quark decays at a future linear e^e~ collider. 

1 Introduction 

One of the most important tasks of current high-energy physics is the search for physics ef- 
fects beyond the Standard Model (SM) , where the Minimal Supersymmetric Standard Model 
(MSSM) [1] is one of the leading candidates. Supersymmetry (SUSY) predicts two scalar 
partners for all SM fermions as well as fermionic partners to all SM bosons. Another impor- 
tant task is the investigation and identification of the mechanism of electroweak symmetry 
breaking. The most frequently investigated models are the Higgs mechanism within the SM 
and within the MSSM. Contrary to the case of the SM, in the MSSM two Higgs doublets 
are required. This results in five physical Higgs bosons instead of the single Higgs boson in 
the SM; three neutral Higgs bosons, /i„ [n = 1,2,3), and two charged Higgs bosons, iJ^. 
In the MSSM with complex parameters (cMSSM) the three neutral Higgs bosons mix [JHS] , 
giving rise to the states hi,h2,h^. 

An interesting production channel of Higgs bosons is the decay of the heavy scalar top 
quark to the lighter scalar top (scalar bottom) quark and a neutral (charged) Higgs boson. 
Another SUSY particle that can produce a Higgs boson is a chargino, which can decay to a 
lighter chargino (a lighter neutralino) and a neutral (charged) Higgs boson. 

The original heavier SUSY particles can be produced at the LHC, or if kinematically 
allowed at an e+e^ collider. At the ILC (or any other future e+e~ collider such as CLIC) a 
precision determination of the properties of the observed particles is expected [6l[7]. Thus, 
if kinematically accessible, Higgs production via scalar top quark or chargino decays could 
offer important information about the Higgs sector of the MSSM. 

In order to yield a sufficient accuracy, one-loop corrections to the various SUSY decay 
modes have to be considered. For the precise evaluation of the branching ratio at least all 
two-body decay modes have to be considered and evaluated at the one-loop level. We review 
the results for the evaluation of these decay widths (and branching ratios) obtained in the 
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MSSM with complex parameters (cMSSM) We will review the numerical results for 



T{i2^iihn) (n= 1,2,3), (1) 
r(i2 ^ hH+) , (2) 
r{x2^X°H-) (j = 1,2,3), (3) 

r(xr ^ x'lH-) , (4) 

where Xfe denotes the neutralinos, Xj'^ the charginos. The total decay width is defined as the 
sum of all the partial two-body decay widths, which have all be evaluated at the one-loop 
level. 

We also concentrate on the decays of ^27X2^1 9'iid do not investigate i^jxti decays. In 
the presence of complex phases this would lead to somewhat different results. Detailed 
references to existing calculations of these decay widths, branching ratios, as well about the 
extraction of complex phases can be found in Refs. [Hlin]- Our results will be implemented 
into the Fortran code FeynHiggs [TIJHT5] . 



2 The complex MSSM and its renormalization 

All the relevant two-body decay channels are evalulated at the one-loop level, including hard 
QED and QCD radiation. This requires the simultaneous renormalization of several sectors 
of the cMSSM, including the colored sector with top and bottom quarks and their scalar 
partners as well as the gluon and the gluino, the Higgs and gauge boson sector with all the 
Higgs bosons as well as the Z and the W boson and the chargino/neutralino sector. Details 
about our notation and especially about the renormalization of the cMSSM can be found in 

Refs. [iiniiiiiin]. 

An important role play contributions of self-energy type of external (on-shell) particles. 
While the real part of such a loop does not contribute to the decay width due to the on- 
shell renormalization, the imaginary part, in product with an imaginary part of a complex 
coupling (such as At or Mi) can give a real contribution to the decay width. These contri- 
butions (in the following called "absorptive contributions" ) have been taken into account in 
the analytical and numerical evaluation. The impact of those contributions will be discussed 
in Sects, d g] 

The Feynman diagrams and corresponding amplitudes contributing to the various decays 
have been obtained with FeynArts |16| . The model file, including the MSSM counterterms, 
is largely based on Ref. |17j . however adjusted to match exactly the renormalization pre- 
scription described in Ref. [8l[9l[T4j[T5] . The further evaluation has been performed with 
FormCalc |18| . As regularization scheme for the UV-divergences we have used constrained 
differential renormalization jl9) . which has been shown to be equivalent to dimensional re- 
duction [20] at the one- loop level [T8|. Thus the employed regularization scheme preserves 
SUSY [2n[22] . All UV-divergences cancel in the final result. (Also all IR-divergences cancel 
in the one-loop result as required.) 



3 Numerical results for scalar top decays 

The numerical examples are shown in two numerical scenarios, SI and S2, where the pa- 
rameters are given in Tab. [TJ The results shown in this section consist of "tree" , which 
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denotes the tree-level value and of "full" , which is the partial decay width including all 
one-loop corrections. We only show the results for the decay widths, since the size of the 
loop corrections to the branching ratios are more parameter dependent. 



Seen. 


tan/3 


Mh± 










At 


Ab 


Ml 


AI2 




SI 


20 


150 


650 


0.4 mr 




200 


800 


400 


200 


300 


350 


S2 


20 


180 


1200 


0.6 nif^ 


0.8 m,-^ 


300 


1800 


1600 


150 


200 


400 



Table 1: MSSM parameters for the initial numerical investigation; all masses are in GeV. 



The production of fa at the ILC(IOOO), i.e. with = 1000 GeV, via e+e" PJ2 wih 
be possible, with all the decay modes ([1]), ([2]) being open. The clean environment of the ILC 
would permit a detailed study of the scalar top decays. For the parameters in Tab. [T]we find 
a{e~^e~ — > tjta) ~ 1-4 fb, i.e. an integrated luminosity of '--^ 1 ab~^ would yield about 1400 fa- 
The ILC environment would result in an accuracy of the relative branching ratio close to 
the statistical uncertainty: a BR of 30% could be determined to ^ 6% for the m^^ values 
in Tab. [TJ Depending on the combination of allowed decay channels a determination of the 
branching ratios at the few per-cent level might be achievable in the high-luminosity running 
of the ILC(IOOO). 

We show the results for the various decay widths as a function of (pAf The other 
parameters are chosen according to Tab. [1] Thus, within SI we have m^^ +"^t2 ~ GeV, 
i.e. the production channel e'^e^ t[i2 is open at the ILC(IOOO). Consequently, the 
accuracy of the prediction of the various partial decay widths and branching ratios should 
be at the same level (or better) as the anticipated ILC precision. 

In Fig. [1] we show r(t2 — iihi) (first), r(i2 — ^1^2) (second), T{i2 hh^) (third) and 
r(t2 -> biH^) (fourth row) as a function of (pAt for the parameters in Tab. [U where the left 
(right) column displays the (relative one-loop correction to the) decay width. While r(t2 — >■ 
tihi) in S2 is of 0(9 GeV), the other decay widths shown arc of 0(1 GeV). The variation 
with (pAt can be seen to be very large, of 0(50%). The size of the one-loop corrections, 
as shown in the right column, is also sizable, of O(±20%) and exhibit a strong variation 
with ifAf The effects of the "absorptive contributions" are clearly visible, especially for 
^2 ~> tihi. Consequently, the full one- loop corrections must be taken into account in a 
reliable complex phase determination from scalar top decays. 

4 Numerical results for chargino decays 

The numerical examples arc evaluated using the parameters given in Tab. [2] We assume the 
scalar quarks to be sufficiently heavy such that they do not contribute to the total decay 
widths of the charginos. We invert the expressions of the chargino masses in order to express 
the parameters fi and A/2 (which are taken to be real) as a function of m ~± and m^± . This 
leaves two choices for the hierarchy of fi and Ar2: 

5> : /i > M2 ixt more higgsino-like) , (5) 
iS< : /i < M2 ix2 more gaugino-like) . (6) 
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Figure 1: Tree- level ("tree") and full one- loop ("full") corrected partial decay widths are 
shown with if At varied. Also shown arc the full one-loop corrected partial decay widths 
including absorptive contributions ("abs"). First row: r(i2 — 5> tihi), second row: r(t2 — > 
ti/i2), third row: T{i2 -> iihs), fourth row: r(t2 -J> biH+). 
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Seen. 


tan/3 


Mh± 


m~± 




Mr 


Mr 


^/ 


S 


20 


160 


600 


350 


300 


310 


400 



Table 2: MSSM parameters for the numerical investigation of chargino decays; all masses 
are in GeV. 



The absolute value of Mi is fixed via the GUT relation (with IM2I = M2) 

|AfiH^tan2 0^M2«iM2 , (7) 

leaving ipMi a-s a free parameter. 

The values of m-± allow xt'>S or xtxi production at the ILC(IOOO) via e+e~ — > 

xfxi2^ with all the subsequent decay modes to a neutralino and a charged Higgs boson, see 
Eqs. (IH). As for the scalar top decays the clean environment of the ILC would permit 
a detailed study of the chargino decays. For the values in Tab. [2] and unpolarized beams 
we find, for 5> (5<), a{e+e~ xfxj) ~ 4(12) fb, and cr(e+e- ^ xfxl) « 55 (80) fb. 
Choosing appropriate polarized beams these cross sections can be enhanced by a factor of 
approximately 2 to 3. An integrated luminosity of ~ 1 ab~^ would yield about 4 — 12 x 10"^ 
xf events and about 55 — 80 x 10'^ Xi'xT events, with appropriate enhancements in the 
case of polarized beams. The ILC environment would result in an accuracy of the relative 
branching ratio close to the statistical uncertainty, see the previous section. Depending on 
the combination of allowed decay channels a determination of the branching ratios at the 
per-cent level might be achievable in the high-luminosity running of the ILC(IOOO). 

The results shown in this section consist of "tree" , which denotes the tree-level value and 
of "full" , which is the partial decay width including all one-loop corrections. Also shown are 
the full one-loop corrected decay widths omitting the absorptive contributions ("full R"). 
We only show the results for the decay widths, since the size of the loop corrections to the 
branching ratios are more parameter dependent. 

In Fig.Hwejhow T{x2 ^ X?i?-) (first), r(x2- X2H-) (second), r(x2- ^ X3H-) 
(third) and r(x;^ — >• XiH~) (fourth row) as a function of (pui for the parameters in Tab. [21 
where the left (right) column displays the (relative one-loop correction to the) decay width. 
The decay widhts are of 0(0.1 GeV) in the case of X2 — >■ XiH~, about five times larger for 
X2 X2 a-nd a factor of ten smaller for the light chargino decay. For the heavy chargino 
decay a strong variation with ipMi can be observed. The size of the one-loop corrections, 
as shown in the right column, is also sizable in the case of the heavy chargino, between 
—4% and -|-12% and show a non-negligible dependence on tpMi- Again the effects of the 
"absorptive contributions" are clearly visible. Also these loop corrections should be taken 
into account in a reliable complex phase determination in the chargino/neutralino sector. 
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Figure 2: Tree-level ("tree") and full one-loop ("full") corrected decay widths are shown 
with (f>Mi varied. Also shown are the full one-loop corrected decay widths omitting the 
absorptive contributions ("full R"). First row: T{x2 — >■ Xi^^")? second row: r(x^ — >■ 
xlH-), third row: T{x^ xlH'), fourth row: T{x^ ^ x°H-). 
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